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Abstract 

The transition between ballistic and diffusive motion poses difficult problems in several fields of physics. In this work we 
show how to calculate the spectra of the correlation functions between fields of arbitrary spatial dependence as seen by particles 
moving through the fields in regions bounded by specularly reflecting walls valid for diffusive and ballistic motion as well as 
the transition region in between for motiion in 2 and 3 dimensions. 

Applications to relaxation in nmr are discussed. 
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I. INTRODUCTION 



The problem of the transition between baUistic and diffusive motion impacts many fields of physics including charge 
transport in semiconductors, electron transport in mesoscopic systems, fluid flow, heat conduction [1 , light transport 
in random media, [2], [3] relaxation and frequency shifts in magnetic resonance of particles moving in non-uniform 
fields [1], [S]. In this work we concentrate on the latter problem and present solutions in finite regions in 1, 2 and 3 
dimensions valid for ballistic and diffusive motions and the crossover between. 

In these problems, the physical effects of interest depend on the spectra of the correlation functions of various field 
components as seen by the moving particles. 
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where () signifies the ensemble average over the gas particles. 

As shown first by [6] (see also [2], [8]) and then applied by [9], [10] and [11] these correlation functions are given by 
the conditional probability p (7^, rl'Fo, 0), the probability that a particle at I^q at t = will be found at a position 
"r^ at the later time r: 

i?i.2 (r) = (Si (<) 52 (t + r)> = J J d\d^roB,il^o)B2{l^)pi^,to + T\l^o,to)p(^a:to) (2) 

wherep (^Oj ^o) is the probability of finding a particle at ~?^o at time io which is taken as constant (= 1/L in 1 dimension). 
In this paper we consider single-speed transport only. 

We treat the motion as a persistent continuous time random walk with isotropic scattering (in 2 and 3 dimensions) 
and uniform velocity with an exponential distribution of times, i, between scatterings, 

^ it) = -e-'/^' (3) 



II. MOTION IN 1 DIMENSION 

Goldstein |12) has shown that in these conditions the conditional probability distribution, considering the motion 
as the continuum limit of a persistent random walk, satisfies the telegrapher's equation in one dimension. As shown 
by Masolivcr ct al [13 the equation can be solved by separation of variables and the conditional probability with 
reflecting boundary conditions at a; = 0, L is given by 
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I X..0) ^ ^ ^ ^c„, COS (=5) ™h g + 1 sinh U e-^ (4) 



with Sn = •\/l — 4a;^r2 and ojn = miv/L with v the velocity of the particles. Note the change in behavior as s„ goes 
from real to imaginary. The behavior of this solution is discussed in some detail in [T3] and [13] . There is a delta 
function peak leaving the source with velocity, v, This peak diminishes in time as particles are scattered and a wake 
of scattered particles builds up behind the peak. The motion can be followed through successive wall reflections. 
Then according to ^ the field correlation function is given by 



i?i,2 (r) = 51 ^1 (fcn) ^2 [kn) <j cosh ( ^ 1 + ^ sinh ["-^]\ e'"^ (5) 



with 



and 
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The Fourier transform of this gives the desired spectrum which determines physical phenomena. 

,2. 



^1,2 (W) = ^ ^1 (fcn) F2 (fc„) 



(8) 



For not too high uj: oj <C ttu/L this result is consistent with one obtained earlier for the diffusion regime of motion 

m, m- 

This result is consistent with calculations of the position- correlation function, Bi = B2 ^ x for various values of 
normalized mean free path, I' = vt^/L presented in [S], see figure 2 in that paper. Notice as the motions approaches 
ballistic, V » 1, the relaxation shows a resonance behavior when a harmonic of the wall collision frequency coincides 
with the Larmor frequency, which has not been noted previously. The resonance peaks are smoothed by averaging 
over the velocity distribution. Experiments to observe this are under way. 

It may be shown that for a perturbation field with a uniform gradient B = Gx the series (6) can be summed in 
closed form: 
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where 



F (w) = 1 - Im 
L 



tan 



{uTc + i) 



Tf) = — characteristic ballistic time 

V 

In the limit ujTc 1 ^ goes over to the diffusion theory result, see [S]. 



(9) 
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III. MOTION IN 2 DIMENSIONS. 



A. Solution in free space 

The obvious extension of the above ideas is to apply them to higher dimensions by writing the Telegrapher's 
equation (TE) with the appropriate form of V^. This was first suggested by Mark Kac, [16]. However the solution of 
the TE cannot represent a conditional probability in two or higher dimensions as it can take on negative values, [17] 
. Morse and Feshbach fl8 discuss the interesting properties of solutions of the TE (as well the wave equation) in 2 
dimensions. This point seems to have been missed in an otherwise interesting and useful work 19J, which nonetheless 
finds the same solution for the spectrum of the conditional probability function as given by [13^. In a remarkable 
paper Masoliver et al |20j have shown that the conditional probability function for a persistent random walk in two 
dimensions, with the time between scattering distributed according to ^ and a uniform distribution of scattering 
angles (s-wave scattering) for particles starting at the origin of coordinates with velocity v, satisfies the 2 D TE with 
an additional source term 



P{r,t) 



2'Kr 



dr 



(11) 



The source moves along with the unscattered particles and represents the scattered particles that are not accounted 
for in the homogeneous TE. The authors give the solution for motion in the infinite domain as: 



P02 {r,t) = e-"-' 



5 [r — vt) 
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where O is the unit step function. The authors also give the Fourier-Laplace transform of the solution 

Po^i^.s^iL.) = I d\ r pir,t)e('^-^-'')dt 



(12) 



(1 -I- iOJTc)'^ + V^T^ 



-,1/2 



(13) 
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The expression (13) and the equivalent one for three dimensions have been rederived by Kolesnik, [21], who gives a 
general treatment applicable to an arbitrary number of dimensions. 



B. Solution in bounded rectangular region. 

We now extend these results to find the spectrum of the field correlation function in a rectangular region bounded 
by reflecting walls. This solution is presented here for the first time. As we desire the spectrum of a field correlation 
function, we will work with the Fourier transform of P02 t) , where we are using the subscript to denote the free 
space solution, 

P02 (^,r) = -±-^ld'Qj duh (^,^) e-^^-^-'^^) (14) 

where P02 = Pm {^,s = iu^ ( 13 ) 

We use the method of images to find the conditional probability function in the presence of the rectangular bound- 
aries at cc = -izLx/^^ y = zLLy/2. The location of the images is sketched in fig.l. Wall reflections are taken into account 
by considering particles coming from the image sources. The desired conditonal probability function for a restrcited 
rectangular domain (with specularly reflecting walls) is given by the superposition of probability coming from the 
original source and all the images. 

We see that the physical region is repeated periodically throughout the plane and there is one image point in each 
cell (the positions of the images are not the same in every cell). Then we can write the solution as 

P(^,t|7^,0) =^Po2(^- A*,t) (15) 

i 

where T^o.i is the location of the i*^ image, (i = denotes the physical source). 

Fig. 2A shows the behavior of P02 C^jt) , (15) for the parameters, v = 1, L — 2, X — 0.2 and t ~ 2.7ih, where 



tfc = L/2?;. Figure 2B show a cross section along the x axis of the function shown in fig.2A. 
Then ^ 

Ri,2 (r) - (Pi (t) B2 (t + T))^ I I cPrcProBi (^0) ^2 (^) ^ P02 - ^0,., r) p (l^o,^, to) (16) 



where the integrals are taken over the physical cell, ^p(^o,i7^o) 



Integrating over (Ptq in the physical cell means that each image point in each cell will cover its entire cell and the 
sum over images and integration over the physical cell can be calculated by taking the function P02 to be the infinite 
domain function integrated over all space if we continue the field periodically as explained in fig.l. This idea was 
introduced by Wayne and Cotts (see also [7]) 



We can rewrite (|16| as 

i?i,2 (r) = / / dxdy \ \ dxodyoB, (T^o) ^2 (^) P02 (^ - ^0, r) (17) 

J^x-t^y J J ~Ly/2 J -oo J -oo 



where Pi (^0) is the periodic extension of Pi (T^o) beyond the physical cell. 

Writing P02 ("r^ — l^o, t) in terms of its Fourier transform (14) we find for the spectrum of the correlation function 
(Jx , ly are integers) 

^1,2 (^) = E (ti.Jy) P02 (-ti.^ly,^) (18) 

with '(ti^,iy = ^) I Ax.iy = L L /-L^/2 /-l"/2 ^■''0^2/06*'^'""'"'^° Pi o) where the integral is taken over one 

period and 

B2 == 1":^^, lX% dxdyB, (^) e-<^-'«-^). 
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FIG. 1. A) Showing the physical region (shaded) and a few of the periodic repetition cells. A source point and its images are 
shown. We consider the conditional probability function as a wave travelling from the source to the obsrvation point, P. It is 
useful to use reciprocity and consider the wave as travelling from P to the source. The solution is then given by considering the 
wave as travelling in an unrestricted domain and taking the sum of all the waves reaching all the image points. In travelling 
from the physical source to the wall at x=L/2 in the physical case with boundaries, the particles see x increasing to the right. 
Trajectories in this region are replaced by trajectories in the region between the image point and the wall so in this region the 
effective value of x=x must be seen to be increasing going from the image to the wall. Thus x has to be taken as periodic with 
period 2L. An arbitrary field depending on position must be treated in the same way, ( as shown in B) which is a plot of an 
arbitrary field varying with x vs. x (solid line)), and hence must be taken as periodic. 



IV. SOLUTION IN 3 DIMENSIONS 
A. Solution in free space 

Masoliver et al [20] (see also [22]) give the form of the Fourier Laplace transform of the conditional probability 
function in three dimensions as 

Po3(^,s) = [ d\ r p{r,t)ei'^-^-'')dt = r^ . .. (19) 

) J Jo (Q.r. -tan-i(i2^)) 

The same result was later obtained by |21j using a different method valid for any number of dimensions. 
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FIG. 2. A) Showing a snapshot of P ( r , r) for parameters given in the text. The white Hnes are the peak of unscattered 
particles while the brightness of the contrast indicates the build up of the wake of scattered particles. B) A cross section 
along the x axis of the function Pq2 r = 2.7) shown in A). The smaller peaks at ±0.7 are the particles that made a coUison 
with the walls a.t x — ±1. The larger peaks represent the crossing of the two "waves" emanating from the image points at 
(±2, ±2) and are larger because they are the sum of the two "waves" . Solid lines are the result of equation ( 15 1 while the dots 
are the results of a Monte-Carlo simulation. 



B. Solution in a rectangular box 



The same reasoning applies here as in 2 dimensions above, the crucial point being that there is one image point per 
cell. 
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The spectrum of the correlation function is then: 

^1.2 M= E Pw.B2{-ti.,i,,i:)h^{-ti.,i,,i.,s^i^) (20) 



where ti.,i,^ = 



1 i'3L^/2 — > ~ 

Pi.,iy,i. = T7 / / da;odyod2oe'^'-'-'^-^°Bi(^o) (21) 

J-L^/2 J-L„/2 J-LJ2 



-L^/2 J-Ly/2 J-L,/2 

(each integration is over one complete period) and 



1 fL^/2 fLy/2 j-L^/2 

^2 (ti.^.ij = 77 / / / dxdydzB2{-^)e-<:^^^^^y^^-^) (22) 

^ J -L^/2 J -Ly/2 J -L^/2 



V. APPLICATIONS 



In the following we show some illustrative examples of the use of the above technique. 



A. Position auto-correlation functions (Uniform gradient field) 



In Fig. 3A) we show the calculation for the spectrum of the x — x auto-correlation function, determining relaxation 
times for a field with uniform gradient, as is usually assumed in studies of relaxation, for the one dimensional case for 
various values of damping using ([9|. 

The results for 2 and 3 dimensions are: 

2 dimensions: 

^ (-) 

3 dimensions: 

(-) 



For large (A 3> L), we see that there is a series of resonances, which will be broadened by averaging over a realistic 
velocity distribution, and non-specular wall reflections. These are Rabi resonances which occur when a harmonic of 
the periodic motion coincides with the Larmor frequency. Similar resonances have been noted in cylindrical geometry 
in another context [21]. For t^ — I the resonances are still slightly visible, while at Tc = 0.1 we have approached the 
diffusion behavior. In 2 dimensions the resonances are less peaked because closed orbits are less probable. In 3D the 
peaks (even in the undamped case) are reduced to steps. 

In the diffusion limit (green curves) we see the usual behavior at low frequencies approaching a constant as w 
(non-adiabatic region). For higher frequencies the spectrum starts to decrease as (l/w^), (adiabatic regime) and at 
still higher frequencies {ujTc ^ 1) it falls as (1/0;**), the super-adiabatic region [22. For moderate (blue) and light 
(red) damping, the l/cj^ behavior disappears and we go directly from a constant to l/w"* behavior. 
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(24) 
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B. Short range interaction between the walls and the spin of the particles 

This problem has been discussed in [lOj using the conditional probability density from diffusion theory. Here in 
figure 4 we present the results valid for the diffusion and ballistic cases and the transition in between calculated using 



the 2 dimensional conditional probability ( 18 ) taking the interaction to be represented by a effective magnetic field: 
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FIG. 3. Showing the spectum of the autocorrelation function of position x for various values of the coUison time, Tc- vs. 
dimensionless frequency curt. A) For one dimension, B )The same for a square in 2 dimensions, C)For a cube in 

three dimensions 
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FIG. 4. Results of the auto-correlation function of the field given in equation (26 1 for the range of force, A = .005, normalized 
to the size of the box and various valus of Tc . 



{L/2+x)/\ _ g-(L/2- 



-x)/Aj 



(26) 



where the range of the interaction is A. 

Here we see the influence of the shape of the field on the spectrum. 

In the medium and underdamped cases the constant behavior goes over, as frequency increases, to a (l/y/uj) 
behavior before reaching the super-adiabatic region (l/oj^) . The extent of the {l/^/uj) behavior will become larger as 
the range of the force decreases. 

In contrast to the spectrum of the position auto-correlation function, in the case of a limited range force the higher 
frequency results are independent of the damping. 



VI. CONCLUSIONS 



We have found expressions for the spectrum of the correlation function of a pair of fields as seen by particles 
executing a continuous time, persistent random walk, with exponential distribution of fiight times in rectangular 
boundaries with specular reflection and single velocity for 2 and 3 dimensions. The results are valid for all values 
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of the damping parameter A/L = Tc/t}j, allowing calculations in the ballistic and diffusive regimes as well as the 
transition region in between. 

These conditional probability spectra allow the calculation of the spectrum of the correlation functions of fields 
with arbitrary position dependence and thus the study of relaxation and frequency shifts in arbitrary fields. The 
results are valid as long as the trajectories are not influenced by the fields. 

It is expected that the results will find applications in other fields such as transport in mesoscopic systems. 
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